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I . INTRODUCTION 



The Air Force has developed an Air Quality Assessment Model ( AQAM) 
to predict the dispersion of pollutants emitted from typical Air Base 
related sources (Ref. 1, 2 and 3) . That model has also been modified to 
better simulate operations from a Naval Air Station (Ref. 4 and 5) . The 
full model consists of three computer programs: Source Inventory, Short 

Term and Long Term Dispersion. In the Short Term Model hourly averaged 
concentrations are predicted over a grid of receptors. Air-base, air- 
craft, and off air-base (environ) sources are considered and are modeled 
as point, area, or line sources as appropriate. 

In the Short Term Model, the jet exhausts from aircraft during taxi 
and takeoff are treated as line sources with no plume rise. Initial 
plume dimensions must be estimated as well as the distance the jet 
"penetrates" the atmosphere before coming essentially to rest relative to 
the ambient wind. Stationary sources such as engine run up test stands 
and test cells are treated as points sources. Likewise, the initial 
plume dimensions must be estimated and a determination whether plume rise 
is a factor must be made. The model then uses these input data to deter- 
mine an initial dispersion coefficient from which a pseudo upwind point 
or line source is located. The sigma's, or dispersion coefficients, are 
based on the stability class of the atmosphere and the downwind distance 
or travel time for an elevated, relatively low velocity source. There is 
a question then, as to the applicability of these dispersion coefficients 
to turbojet exhausts which are often sonic, have high shear stress/ turbu- 
lent mixing characteristics and have different turbulence scales than 
typical exhaust stacks. Evidence that there is a discrepancy has been 
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indicated by an Air Force study (Ref. 6) which presented a pictorial in- 
vestigation of the effects of atmospheric stability on a jet exhaust. 

It is evident that qualitative and quantitative data are needed to char- 
acterize the dispersion of a turbojet exhaust during ground operations. 
This investigation attempted to determine turbojet exhaust dispersion 
as a function of jet characteristics and wind direction under laboratory 
simulated, neutrally stable atmospheric conditions. 

The first part of the study involved the correct simulation of the 
lower portion of the atmoshpere in an aeronautical type wind tunnel. 
Published studies (Refs. 7, 8 and 9) have indicated that it is not 
sufficient to simulate only the mean velocity profile of the atmosphere. 
The turbulent structure of the atmosphere must also be correctly modeled 
simultaneously with the mean velocity profile. Proper simulation can be 
accomplished following the work of Sundaram and Ludwig (Ref. 7) , tripping 
the boundary layer with a suitable fence at the entrance of the wind 
tunnel test section and then letting the turbulent flow develop over a 
length of roughness elements. Suitable matching between the boundary 
layer trip and the roughness elements then results in an area of constant 
shear stress which adequately models the atmospheric flow. 

Once an atmoshperic surface layer has been suitably simulated, a 
turbojet exhaust must then be simulated. If it is assumed that heat and 
mass have the same turbulent diffusion characteristics, then the disper- 
sion rates of pollutants can be determined by measuring the temperature 
distributions in the vertical and horizontal directions. 

In this investigation data were obtained for varying angles of inci- 
dence to the oncoming wind for choked and unchoked nozzle flow conditions. 
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for two jet temperatures r and for two conditions of nozzle surface block- 
age. The data were compared to the dispersion routines used in the AQAM 
computer program. 
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II. ATMOSPHERIC SIMULATION 



In order to make meaningful wind tunnel measurements to determine 
jet aircraft exhaust dispersion rates during taxi, idle and takeoff 
modes of operation it was necessary to correctly model the lower struc- 
ture of the atmoshpere. Sundaram and Ludwig (Ref. 7) have published a 
theoretical study of the requirements for the modeled flow and an ex- 
perimental study of the flow generated by wind tunnel techniques. This 
reference was used as the principal guide in modeling the atmoshperic 
surface layer in a wind tunnel at the Naval Postgraduate School. Sup- 
plementary discussions are presented in Refs. 10-12. 

A. THEORETICAL REQUIREMENTS FOR ATMOSHPERIC MODELING 

Reference 13 describes the atmoshperic boundary layer as a turbulent 
layer which is influenced by a combination of actions due to surface 
friction, thermal stratification, and Coriolis force. The atmospheric 
boundary layer is divided into two distinct layers. The lower region 
called the surface boundary layer extends up to a nominal height of 100 
meters. The region above the surface layer is called the Ekman layer 
and extends up to a height of about 1000 meters. In the surface layer, 
the wind direction is generally considered to remain constant and the 
Coriolis force is not an important parameter. In an ideal environment 
the near-surface region is horizontally homogeneous and free of pressure 
gradients. These conditions imply that the shear stress can be considered 
to be independent of height in the region close to the surface. Surface 
layer theories (Ref. 14) often employ the assumption that the flow char- 
acteristics close to the ground can be expressed entirely in terms of 
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conditions at the surface. Most atmoshperic boundary layers are observed 

to have no laminar sublayer near the surface and are described to be fully 

aerodynamically rough. A turbulent, neutrally stable flow structure in 

the atmospheric surface boundary layer can be described entirely by the 

following parameters; the kinematic viscosity v , the friction velocity 

u. , and the roughness parameter z (Ref. 7) . u. is defined as being 

o 

equal to ( t /p ) and z is a length-scale which describes the influ- 
o o 

ence of a rough surface on a particuJ ar flow. 

Using the above parameters, an aerodynamically rough flow exists if 
(Ref. 7) 

z u /v>3. (1) 

o * 

If z^ is << than the vertical dimension of the boundary layer, then the 
mean- velocity profile of the wind is given in logarithmic form by the re- 
lationship (Ref. 7) 

u/u* = (1 /k) In (z/z ) , (2) 

o 

where u is the average velocity, z is the vertical distance from the 
surface and k is Von Kantian' s constant. This equation is a form of the 
well known law of the wall and is indicative of the interaction between 
the air flow above a given surface and that surface. 

To accurately model the atmosphere the requirements of aerodynami- 
cally rough flow, horizontal homogeneity and the absence of pressure 
gradients must be adhered to. Strict horizontal homogeneity, however, is 
extremely difficult to obtain in the laboratory. Partial simulation, 
therefore, generally has to be imposed. A complete discription of the 
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effects of partial simulation is given in reference 7 and will not be 
repeated here. The conclusion reached from reference 7 regarding partial 
simulation was that a laboratory flow which was fully aerodynamically rough, 
nearly horizontally homogeneous and relatively free from pressure 
gradients was a reasonable model for a neutrally stable atmosphere. The 
generation of such a flow is described in the next section of this report. 

B. EXPERIMENTAL REQUIREMENTS TO SIMULATE ATMOSPHERIC FLOWS 

Aircraft jet exhaust dispersion in the "near ground" environment was 
the primary concern of this investigation; therefore only the surface lay- 
er of the atmosphere was required to be simulated. As discussed previous- 
ly, this was possible since both the atmospheric flow and the simulated 
flow are described completely in terras of conditions at the surface (ie. 

the parameters u. , z , and v ) . The similarity requirements were, 

* o 

therefore, not a function of the geostrophic wind in the real flow nor a 
function of the free stream velocity in the simulated flow. Horizontal 
homogeneity cannot be strictly satisfied since for an experimental flow 
over a flat plate, equilibrium conditions change continually as a function 
of downstream distance (Ref. 7) . 

Since the atmosphere is assumed to exhibit a logarithmic mean velocity 
profile and to be horizontally homogeneous, this part of the similarity 
requirement can be verified by measuring the mean velocity profile at 
several locations in the wind tunnel flow. In reference 7 the turbulence 
intensity and the integral scale of turbulence were also measured and com- 
pared to available atmospheric data. Turbulence intensity measurements 
were also made in this experiment and compared to the results of reference 
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7. Integral scales of turbulence were not measured. For virtually 
identical experimental set-ups it was assumed that if all other experimental 
measurements compared favorable to those of reference 7, then the scales 
of turbulence would also compare favorably. 

C. EXPERIMENTAL APPARATUS AND MEASUREMENTS 

1. Wind Tunnel Modifications 

The experiment was conducted in a low speed wind tunnel at the 
Naval Postgraduate School (NPS) , Monterey, CA. The wind tunnel draws am- 
bient air from the surrounding area into a 5X5 foot test section approx- 
imately 21 feet long. In order to be comparable with the wind tunnel used 
in reference 7, the wind tunnel at the Naval Postgraduate School was modi- 
fied in the following manner. A false ceiling was installed to maintain 
an approximate zero pressure gradient in the axial direction. A two inch 
high boundary layer trip constructed of wood was placed across the entrance 
to the test section. The trip was followed by 2.5 X 2.5 feet square 
sections of roughness elements which covered the entire floor of the test 
section. The roughness elements were 4 sided regular pyramids which were 
0.75 inches square at the base and 0.75 inches in height. They were con- 
structed out of 1.5 inch thick fiberboard. The tunnel and roughness ele- 
ments are shown in figure 1. 

2. Pressure Measurements 

Six, 1/8 inch diameter pressure taps spaced 4 feet apart, were 
mounted flush along the centerline of the wind tunnel wall. The differen- 
tial area required to maintain a zero pressure gradient in the axial direc- 
tion was determined from an initial set of pressure readings. The false 
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ceiling which implemented the required area change was then installed and 
the pressure readings were again taken to verify the required zero axial 
pressure variation. 

3. Electronic Equipment 

All velocity related measurements were taken with a Thermo Systems 
Inc. linearized hotwire anemometer system. DISA single wire and cross wire 
probes and probe holders were used with the anemometer system (fig. 2 and 3) . 
A DC digital voltmeter and a true RMS meter were used to record the steady 
and fluctuating single wire probe data respectfully. A model 1015C Thermo 
Systems Inc. correlator was used in conjunction with the cross wire outputs 
to obtain the sum and difference of the two signals (fig. 4) . The true 
RMS meter was also used to record the cross wire outputs via the correlator. 
Tne hotwire probes were mounted under the wind tunnel such that the probes 
traversed the wind tunnel in a vertical direction. A Spectral Dynamics, 

Spec tra-S cope Model SD-330 frequency counter was used in conjunction with 
the measurement of the centerline velocity (fig. 5) . 

4 . Wind Tunnel Centerline Velocity Measurements 

The maximum velocity expected in the wind tunnel was between 40 
and 50 feet per second. Existing pitot-static systems at NPS were not sen- 
sitive to this small velocity. Therefore, the relationship between the 
Strouhal Number and the Reynolds Number was used to determine the center- 
line velocity (Ref. 15) . The Strouhal number relates the velocity of a 
flow to the frequency at which vortices are shed behind a circular cylin- 
der which is immersed in that flow. 

Sn (Strouhal Number) = (freq. X diam. of cyl.) / velocity. 
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A 0,075 inch diameter stainless steel wire was attached to a collar which 
was fitted over a single hotwire probe such that the hotwire was in only 
one wake (shed vortex sheet) of the cylinder (fig. 6) . The probe, collar 
and hotwire combination was placed perpendicular to the flow at the center- 
line of the tunnel. With the tunnel on, the hotwire system was sensitive tc 
the velocity of the free stream and the output was observed and recorded 
on the frequency counter. The observed frequency of shed vorticies and 
the unique relationship between the ^trouhal Number and the Reynolds Num- 
ber were used to iteratively determine the velocity of the free stream. 

5. Velocity Profile Measurements 

The mean velocity profile measurements were made using the single 
hotwire apparatus. The hotwire system was calibrated in the linear mode 
such that the centerline velocity was made equal to one volt. The single 
wire was also calibrated for flow direction sensitivity. This was 
accomplished by rotating the probe until a maximum reading was observed on 
the DC digital voltmeter. It was then assumed that the probe was perpen- 
dicular to the flow. All subsequent measurements were referenced to this 
direction. The linearized output was averaged through a ten second time 
constant before being displayed on the DC digital voltmeter. The velocity 
profile was obtained by recording the DC value of the hotwire output as it 
traversed the wind tunnel in the vertical direction. Vertical measurements 
were referenced to the top of the roughness elements and were measured to 
within an accuracy of 0.020 inches with respect to that reference. Veloci- 
ty profiles were obtained at two axial locations. The locations were 72 
inches and 109 inches downwind from the boundary layer trip at the entrance 
to the test section. The two axial measurements were necessary to ensure 
the required degree of uniformity in the jet dispersion region. 



9 



6. Turbulence/Shear Stress Measurements 



The cross wire system described earlier was used to measure the turbu- 
lence intensities and shear stresses directly. The cross wire was calibra- 
ted such that the linearized output for each wire was of equal sensitivity. 
This was accomplished by first orienting one wire normal to the flow and 
then adjusting its sensitivity to the proper level. The probe was then 
re-oriented so that the second wire was normal to the flow and its sensi- 
tivity was then matched to that of the first wire. For data acquisition, 
the probe was oriented axially such that the longitudinal turbulence in- 
tensity was proportional to the RMS value of the sum of the cross wire 
outputs. Likewise, for the same orientation, the vertical turbulence in- 
tensity was proportional to the RMS value of the difference of the cross 
wire outputs. The turbulent shear stress is defined as - pu^w' . This 
shear stress was obtained directly since u'w'' is proportional to the 
square of the sum minus the square of the difference of the cross wire 
outputs. Like the single wire data, the cross wire data was recorded as 
the probe was traversed in the vertical direction at the same two axial 
locations. For redundancy, the vertical profile of turbulence intesnity 
was also measured using the single wire probe. 

D. RESULTS AND DISCUSSION 

All the measurements obtained at the two axial locations, ie. 72 and 
109 inches from the entrance of the wind tunnel test section, were taken 
with a wind tunnel center line velocity of approximately 40 feet per 
second. 
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1 . 



u x and z 
o 

From the law of the wall, 

u/u. =5.75 log z/z (3) 

* o 

Equation (3) can be rewritten as 

u/U = (u./U ) 5.75 (log z - log z ) . (4) 

00 " 00 Q 

When u/U^ is plotted vs. log z, 

slope = 5.75 (u*/Uj . (5) 

u/U m vs. log z was plotted for both axial positions and a least 
squares fit was obtained through the data points (figures 7 and 8) . The 
slopes were measured and u* was determined by; 

u. = (slope • U )/5.75. (6) 

u* values of 2.32 and 2.52 feet per second were obtained at the 72 and 
109 inch positions respectively. Using the least squares fit data, the 
calculated u. was used to determine z by; 

log z = log z - u/(5.75 • u.) (7) 

o * 

z values of 0.0215 and 0.0326 were obtained at 72 and 109 inches 
o 

respectively. 

2 . Vertical Distribution of Shear Stress 

Since u* is proportional to the shear stress, 

u. = /t /p = - \/u' w' , (8) 

* o 

a measure of the average shear stress for each position can be obtained 

2 

from the mean velocity profile data. Non-dimensionalizing by ; 

u ' w ' / U 2 = (u */ U M )2 
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The measured shear stress, u'vr' , from the cross wire data, was also non- 

2 

dimen tionali zed by U and plotted as a function of the vertical dis- 

oo 

tance above the roughness elements for both axial positions. Figure 9 
presents the cross-wire data. Also shown are the values of average stress 
determined from the law of the wall. 

Figures 10 and 11 compare the results with those presented in refer- 
ence 7. It is observed that the shear stress continued to increase until 
very near the wall whereas the shear stress from reference 7 began to de- 
cay further from the wall. Figure 10 shows that the shear stress was 
approximately constant for heights from 2 to 5 inches above the rough- 
ness elements at the 72 inch position. At the 109 inch position 
(fig. 11) , the shear stress decayed gradually from the near wall region. 

In an attempt to modify the near wall shear stress distribution, different 
combinations of trip heights and wind speeds were tried. This proved to 
be unsuccessful. The mean velocity profile was observed to vary much more 
rapidly in the axial direction and the relatively constant shear stress 
region obtained with the two inch trip was destroyed in varying degrees 
depending on which trip or what tunnel speed was used. There was also a 
somewhat greater variation in horizontal homogeneity in this experiment 
relative to reference 7. This was reflected in the increased slope of 
figure 8 over that of figure 7. This was also reflected in figure 9 by 
noting that the average shear stress differed by about 2 0 percent 

between the two positions. Moreover, z increased by about 50 

o 

percent over the axial range compared to about 30 percent for the experi- 
ment of reference 7. It was also noted in reference 7, however, that the 
degree of uncertainty in the atmosphere may be much larger than the 30 
percent figure. 
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3. Vertical and Horizontal Turbulence Intensities 



Figures 12 and 13 compare the measured vertical and honzont* 
turbulence intensity measurements at the two axial locations with 
of reference 7. In general, the vertical and horizontal data ag 
tremely well with the data of reference 7. The horizontal turbul 
a little lower for both positions in the range of Z / Z Q from 1 
40. This corresponded to a a height of about one inch above the 
elements and was the area of difference noted in figures 10 and 11 
Since, (1) a region of constant shear stress was obtained i 
agreement with reference 7, (2) the turbulence intensities were 
constant over the same region, (3) the roughness parameter z 

o 

parable to that of reference 7, and (4) equation (1) was satisfied, 
simulated atmoshperic surface layer was considered to be adequate 
initial study of jet exhaust dispersion. 
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III. SIMULATED JET DISPERSION IN A WIND TUNNEL 



A. BACKGROUND 

The AQAM model treats most aircraft emissions as finite line source 
Aircraft engine test cells and run up stands are treated as point source < . 
Line sources are modeled by a general line source or puff type model 
(Ref. 1-3) . Point sources are modeled by the application of Gaussi^ 
plume theory. This theory represents the downwind concentration of pol 
lutants from an elevated source as Gaussian distributions in both the 
horizontal and vertical directions. When applied to aircraft exhaust: , 
the line source method first assumes that the emissions from the moving 
aircraft come to rest relative to the ambient wind after a penetration 
length due to the velocity effects of the jet. After coming to rest, the 
source of finite cross-section is generally segmented. The segmented 
portions of the source are then assumed to behave as pseudo-upwind line 
sources which are dispersed downwind in a Gaussian manner. Point and 
line sources both employ the same dispersion coefficients a ^ and 
which are the standard deviation points of a Gaussian or normal distn 1 
tion curve and are described by Turner in his workbook on dispersion 
estimations (Ref. 16) . As mentioned earlier, the dispersion coeffici "t 
are based on the atmospheric stability class and the time or distance 
traveled downwind from a relatively low velocity source. The applicabilii 
of these coefficients relative to a high velocity, horizontally emitted 
source, such as a turbojet engine, was investigated for various jet para- 
meters. 
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B. EQUIPMENT 



1. General Discussion 

The investigation was conducted in the low speed wind tunnel 
described in the previous section, A bumer/nozzle assembly provided the 
simulated turbojet exhaust. The temperatures within the exhaust were 
measured in the vertical and horizontal directions using a thermocouple 
apparatus . 

2 . Burner/Nozzle 

Oxygen and ethylene were mixed and ignited at the entrance to a 
burner can, figure 14. The hot gases were diluted with air to provide 
burner cooling and the desired flow rates. The burner exhausted into the 
wind tunnel through a 1/8 inch diameter stainless steel tube. The tube 
was inserted through the tunnel floor. It was made with one 90 degree 
bend which aligned the exhaust parallel to the tunnel floor. A nozzle was 
created by reducing the exit diameter of the tubing to 0.058 inches. The 
tube was positioned 3 inches above the roughness elements. This position 
was chosen since it was about the center of the area of constant shear 
stress discussed in the atmospheric simulation section of this report. 
Stagnation pressure and temperature probes were used to determine the 
nozzle exit conditions. The nozzle entered the wind tunnel at a position 
72 inches downstream from the boundary layer trip described in the previ- 
ous section. 

3 . Temperature Measuring 

Horizontal and vertical temperature profiles were obtained using a 
chromel-alumel thermocouple probe. The thermocouple was referenced to an 
Omega electronic Miniature Cold Junction (MCJ) . The nozzle and 
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thermocouple system is shown in figure 15. A probe holder and a traverse 
system permitted positioning of the thermocouple probe in varying axial, 
transverse, and vertical locations. The traverse system, an electrical- 
mechanical device mounted under the wind tunnel (fig. 16) , facilitated 
positioning in the vertical and transverse direction. The probe holder, 
mounted on the traverse system, permitted the thermocouple probe to be 
positioned in the axial direction. Two X-Y plotters (fig. 17) were 
used to record the temperature profiles as a function of downwind distance. 
One plotter was used for the vertical temperature profile and the other 
for the horizontal temperature profile. 

C. DISCUSSION OF RESULTS 

1. Introduction 

The dispersion data from the experimental jet exhaust were used to 
determine 1) if the plume dispersion was Gaussian, 2) the sensitivity 
of the jet dispersion to the wind direction, 3) the sensitivity of the 
jet exhaust dispersion to nozzle exit Mach number (choked and unchoked) 
under no-wind and wind conditions, 4) the sensitivity of the jet disper- 
sion to upwind nozzle surface blockage and 5) the effect of increased 
jet exhaust exit temperature on plume rise. The first four items were in- 
vestigated with an exhaust stagnation temperature of 550 degrees F while 
the stagnation temperature for the plume rise investigation was increased 
to 700 degrees F. For choked flow the exhaust stagnation pressure was 32 
psia. For unchoked flow it was 22 psia. 

2 . Attainment of Gaussian Dispersion 

A Gaussian or normal distribution with a mean value of zero is 
defined by the equation 
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y = {1/ ( /2 ttc) } exp [-.5 (x/a) 2 ] 



( 10 ) 



A standard normal distribution (fig 0 18) is a normal distribution in which 
a is set equal to one. In this case, the maximum value of the ordinate 
of equation 10 is l//2rr or 0.3989. The value of the ordinate at one 
standard deviation, where x = c , of equation 10 is 0.242. To determine 
whether the measured temperature profiles were Gaussian the profiles were 
standardized in the following manner. For any measured temperature profile 
the peak temperature (Tmax) was assumed to be at the mean of the distri- 
bution. That temperature was scaled such that Tmax/C was equal to 0.3989. 
The ordinate values of the temperature profile were then scaled by the 
factor C. The ordinate value of the one a point, T^ , of the profile 
was determined by; T^ = C# 0.242. At this temperature the width of the 
temperature plume is equal to 2a . The standard deviation which was 
determined in this manner was used to non-dimensionalize or standardize 
the abscissa of the temperature profile. The standard deviations are ta- 
bulated in Tables I, II and III for each of the conditions of the experi- 
ment. The standardized data points were then plotted against a standard 
normal distribution curve and compared for closeness of fit. Figures 19 
through 24 compare the vertical and horizontal temperature profiles for a 
choked nozzle flow to the standard normal distribution curve. The jet 
exhaust was oriented at zero degrees incidence to the free stream wind 
direction. The wind tunnel centerline velocity was 40 feet per second. 

The different sets of data represent measurements taken at different axial 
positions downwind from the jet exhaust. The downwind distance was ex- 
pressed in terms of jet diameters. The plume was near Gaussian after 
approximately eight jet diameters downstream from the exit plane of the 
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jet for both the horizontal and vertical profiles. The distance required 
for the profile to become Gaussian appeared to be the result of the expan- 
sion process of the fluid as it exited the nozzle. Similar plots were ob- 
tained for the other zero incidence conditions shown in Tables I and II. 

For each case the plume exhibited near Gaussian properties at about eight 
jet diameters downstream. The dashed line on Tables I and II indicate the 
position at which the plume exhibited Gaussian characteristics. 

As the nozzle was rotated to the wind direction, the jet-wind inter- 
actions caused the horizontal profile to be negatively skewed. The degree 
of skewness diminished as a function of downwind distance as the effect of 
the jet exit velocity was overcome by the force of the wind. At the point 
where the wind was the dominant dispersion factor (ie. past the penetration 
length) the plume again exhibited Gaussian characteristics. Figure 25 com- 
pares the data obtained at an axial distance of 1.72 jet diameters for 
different angles of incidence to the wind. The angles investigated were 
29, 45, 67 and 90 degrees and the data were non-dimensionalized in the 
same manner as discussed previously. The degree of skewness in figure 25 
is evident for each angle of incidence when data are compared to the stan- 
dard normal distribution curve. Similar data for different downwind sta- 
tions are plotted on figures 26 through 32. The vertical temperature 
profiles exhibted no skewness and were very nearly Gaussian at each data 
position. The a's which were obtained for the vertical and horizontal 
dispersions are tabulated in Table III. The downwind distance at which 
the plume became Gaussian is indicated for each angle by a dotted line. 
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3 . Sensitivity of Exhaust Dispersion to Angle 



of Incidence to the Wind 

A spatial representation of the jet exhaust as a function of du 

i 

wind and horizontal distances from the centerline of the nozzle is plot 
ted on figure 33 for jet exhaust angles of incidence to the wind of 29, 
45, 67 and 90 degrees. As expected, as the angle of incidence was in 
creased, the length normal to the wind direction where the dispersion 
tribution was essentially uniform (ie. sigma approaches infinity) incr 
(in a nonlinear relationship, fig. 34) . The downwind distance at whir 
the plume became approximately uniformly distributed was extremely sens 
tive to the angle of incidence. The downwind distance decreased quite 
rapidly as the angle of incidence was increased in an almost linear re- 
lationship (fig. 34) . From Table III it can also be seen that the 
corresponding sigmas increased considerably. Quantitative data was 
attempted for angles of incidence greater than 90 degrees. The probe 
could not be translated further upstream than the 90 degree jet orient 
tion. At this position, the exhaust of the jet was completely disper. 
Qualitatively then, as the angles of incidence were increased, the exit 
jet momentum and wind action acted together to rapidly disperse the ex- 
haust in the downwind direction. The dispersion coefficients were in- 
creased which resulted in a decreased downwind distance where the plume 
became uniformly distributed. 

4 . Sensitivity of Dispersion Rate to Nozzle Exit Mach Number 
(choked and unchoked) , for Wind and No-wind Conditions 
Table I presents the results in terms of the size of the one 
standard deviation point of the plume in jet diameters as a function of 
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downwind distance from the nozzle exit. It can be seen that the plume di- 
mensions were relatively insensitive to varying exit conditions (choked, 
unchoked and high temperature) . The data for the unchoked nozzle-with-wind 
dispersed about the same as did the choked flow. For choked flow the 
Mach number at the exit was 1.0. For the unchoked flow the Mach number 
was about 0.8. Therefore, both had high velocities relative to a typical 
stack exhaust. 

The dispersions from a choked nozzle with and without wind were also 
measured. The two sets of data for both the horizontal and vertical di- 
rections are plotted in figures 35 and 36. It is seen that the disper- 
sion of the plume became affected by the wind structure at about 30 jet 
diameters downwind from the nozzle exit 0 Between the nozzle exit and 
about 30 jet diameters downstream, the no-wind and wind data compare 
fairly well with each other. Further downwind from that position, however, 
the two curves diverge. It was apparent then that the exit conditions of 
the nozzle (ie. the penetration length) governed the dispersion of the 
plume for about the first 30 jet diameters downwind, and thereafter the 
dispersion was primarily a function of the environmental conditions. 

5. Effect of Nozzle Surface Blockage 

To investigate the effect of flow disturbances (caused by the 
tubing upwind of the nozzle) on the dispersion of the jet exhaust, two 
different non-heat conducting disks were separately fitted on the nozzle 
tubing. These blockage devices increased the upstream area relative to 
the 1/8 inch tubing by a factor of 15 and 75 respectfully. Table II 
contains the results and it was observed that upstream blockage did not 
greatly affect the dispersion rates. These upstream blockage effects 
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were apparently dissipated by the time that the jet had penetrated approx- 
imately 30 jet diameters aft of the nozzle. These data indicate that air- 
craft configuration may not significantly affect the jet dispersion for 
the slowly moving taxi operation. 

6. Sensitivity to Increased Exit Temperature 

As stated previously, the effect of exit temperature had little 
effect on the plume dimensions as a function of downwind distance. An 
investigation of the plume rise was conducted by measuring the vertical 
distance between the peak point of the temperature profile and the center- 
line of the nozzle for two exhaust temperatures . The two sets of data are 
compared as a function of downwind distance in figure 37. Because of the 
small distances it was difficult to make accurate quantitative conclusions 
with respect to the plume rise. From the data, however, it appeared that 
the plume rise was about one jet diameter per 100 degrees F change in exit 
temperature. Equally as interesting as the plume rise was what appeared 
to be a small degree of looping caused by the turbulent structure of the 
atmosphere. This phenomenon is shown in figure 37 as a dotted line. 

D. COMPARISON OF EXPERIMENTAL RESULTS WITH AQAM 

The line source dispersion employed for aircraft exhaust jets in AQAM 
first assumes that the jet exhaust has a penetration length equal to 140 
meters (ie. approximately 140 jet diameters) . The penetration length is 
that length where the velocity effects of the exhaust cease to be a fac- 
tor in dispersing the exhaust. Thus, the exhaust essentially comes to rest 
with respect to the ambient air mass. At this distance then, the exhaust 
is being dispersed entirely by the action of the wind. 

In AQAM, estimates of the lateral and vertical dimensions of the 
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aircraft plume are made at the penetration length. The lateral width w , 

and vertical height Az then define the dimensions of a line source of 

length L . The lateral width is assumed to remain constant over the 

length of the line. Once the line source has been defined, the line is 

generally segmented and dispersed downwind using the Gaussian theory of 

dispersion. That is; 'from the dimensions of the line source an initial 

o and a are determined from empirical relationships 0 = w/2.4 and 

yo zo yo 

a = Az/2.4 respectively. With these values of a and a an up- 
zo yo zo 

wind distance for a pseudo upwind line source is determined using the 

a and a vs. distance or time charts. The pollution concentrations 
yo zo 

are then dispersed from the determined line source in a Gaussian manner. 

If w and Az are not specified inputs, the AQAM model uses default 
values of 20 meters and 8 meters respectfully. 

The penetration length determined in this investigation for zero 
degrees orientation to the wind was approximately 30 jet diameters, con- 
siderably shorter than the 140 jet diameters used in AQAM. The 30 jet 
diameter figure was further investigated using the dispersion data from 
tests with the nozzle oriented at other than zero degrees to the wind 
(fig. 38) . The 30 jet diameter penetration distance determined from the 
axial conditions, indicated by an arc in figure 38, was compared to the 
point where the dispersion curves exhibited no skewness (ie. where the 
distributions were Gaussian normal to the wind direction) . It was 
assumed that beyond this position the jet effects could be neglected and 
that the jet would be disperesed by the ambient wind only. This 
assumption is not entirely correct since the plume continued to exhibit 
horizontal displacements normal to the wind. Along the jet axis the 
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penetration length increased as a function of incidence to the wind from 
about 28 to 58 jet diameters, still considerably less than the 140 jet di- 
ameters assumed by AQAM. This relationship is plotted in figure 39 and 
appears to be approximately linear between 29 and 90 degrees relative to 
the axial orientation. The distance downwind from the nozzle where the 
plume became Gaussian, however, remained approximately constant at 25-30 
jet diameters. It was also noted that the corresponding a's at the 
Gaussian point for each angle became wider with increased angles of inci- 
dence due to more turbulent mixing. This relationship is plotted in 
figure 40. With the limited number of data points available it appears 
that the relationship follows a smooth curve. a ^ and and the jet 

penetration length seem to be predictable between zero and 90 degrees of 
incidence to the wind for neutral stability and these nozzle exit con- 
ditions. These values could be used in AQAM to determine w and Az 
instead of estimating them or using the default values. The relationship 
shown in figure 39 could also be used as input for jet penetration length. 
However, since the penetration lengths are small, it may be better to 
neglect this effect entirely. 

A comparison of the rates of dispersion for zero degrees incidence 

to the wind (i.e. the slopes of the a and a vs. distance curves) was 

y z 

made between those for a neutrally stable atmoshpere (class D) in 
reference 1 and the experimental data. As can be seen from figures 35 
and 36 the experimental dispersion rates were greater than the class D 
rate used by reference 1 for a neutrally stable atmosphere. Reference 1 
indicated that this may be expected. The dispersion rates were obtained 
from data compiled from typical low velocity stack emissions and not from 
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high velocity sources. In general terms, it appears that for the condi- 
tions of this experiment, a more realistic representation of the actual 
dispersion rate could be obtained if the atmospheric stability class was 
decreased by one. 

For a vertical stack, AQAM also estimates the plume rise due to 
thermal effects by one of two relationships. These relationships were 
investigated to see if they could be used to predict the plume rise ob- 
served in the experiment. Both Holland and Carson-Moses formulas 
were found to predict insignificant plume rise for the conditions of this 
experiment. 
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IV. SUMMARY AND RECOMMENDATIONS 



in summary, a neutrally stable atmospheric surface layer was simu- 
lated in a low speed wind tunnel. Horizontal and vertical temperature 
profiles aft of a converging nozzle were measured. The dispersion rates 
from the simulated jet exhaust were measured for different nozzle 
conditions and compared to the dispersion rates used in the AQAM Model. 

Shorter jet penetration lengths and more rapid jet dispersions were 
observed. Initial plume dimensions were found to vary significantly with 
jet orientation to the ambient wind direction and some plume rise was 
observed. 

It is recommended that additional studies be made for conditions of 
other than neutral stability. Measurements made in less stable surface 
layers would provide a better understanding of the plume rise and looping 
phenomena. 
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Downwind Distance JD 





Sigma Vert, 
(choked) 
Blockage 


(JD) 


Sigma Horz. 
(choked) 
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3.01 


34.5 


2.90 


3.16 


3.44 


3.01 
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3.87 


51.7 
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4.59 


4.59 


4.74 
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7.75 


7.32 


7.50 


7.75 


120.6 


9.19 


9.19 


9.19 


11.20 


11.30 


10.34 



Table II 



1-Sigma Plume Width of the Vertical and Horizontal 
Temperature Profiles at Zero Degrees 
Incidence to the Wind for Different 
Nozzle Blockages 
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1-Sigma Plume Width of the Vertical and Horizontal 
Temperature Profiles for Different Angles of 
Incidence to the Wind 
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Figure l 
Single Hotwire 







Figure 3 
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Figure 4 

Anemometer and Correlator 




Figure 5 

Frequency Counter 





Figure 6 
Velocime ter 
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Figure 7 

Velocity Profile at 72 inches Downwind 
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Figure 8 

Velocity Profile at 109 inches Downwind 
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Figure 9 

Vertical Distribution of Shear Stress at 
72 and 109 inches Downwind 
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Height above Roughness Elements, Z (inches) 
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Figure 10 

Vertical Distribution of Shear Stres* at 
72 inches Downwind Compared to the Data of Ref. 7 
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Figure 11 

Vertical Distribution of Shear Stress at 
109 inches Downwind Compared to the Data of Ref. 7 
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Comparison of Turbulence Intensities at 
72 inches Downwind to the Data of Ref. 7 
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Figure 13 

Comparison of Turbulence Intensities at 
109 inches to the Data of Ref. 7 




Figure 14 
Burner Can 
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Figure 15 

Exhaust Nozzle and Thermocouple Probe 




Figure 16 

Thermocouple Traverse System 
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X-Y Plotters 



Figure 17 

for Temperature Measuring 
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Standard Normal Distribution 



Downwind Distance 
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Comparison of Vertical Distribution of 
Temperature to a Standard Normal Distribution 
for Zero Degrees Angle of Incidence 
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Comparison of Vertical Distribution of 
Temperature to a Standard Normal Distribution 
for Zero Degrees Angle of Incidence 
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Comparison of Vertical Distribution of 
Temperature to a Standard Normal Distribution 
for Zero Degrees Angle of Incidence 
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Downwind Distance 
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Figure 23 

Comparison of Horizontal Distribution of 
Temperature to a Standard Normal Distribution 
for Zero Angle of Incidence 
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Angles of Incidence 

A- 29 Degrees 
• - 45 Degrees 
■ - 67 Degrees 
V - 90 Degrees 
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Comparison of Horizontal Distribution of Temperature to 
Standard Normal Distribution for Angles of Incidence at 
1.72 Jet Diameters Downwind from the Nozzle 



Angles of Incidence 

A - 29 Degrees 
• - 45 Degrees 
■ - 67 Degrees 
V- 90 Degrees 




c0 


4J 




o 


CO 


a> 


4-J 




rH 




a) 


N) 


Q) 


o 


N 


U 


c 


O 


3 


a> 


s 


4J 







◄ 



i H 

ro 

6 o 



cvj — : O 

6 o 



O 



52 



Comparison of Horizontal Distribution of Tempera 
Standard Normal Distribution for Angles of Incii 
8.62 Jet Diameters Downwind from the 
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Angles of Incidence 
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Comparison of Horizontal Distribution of Temperature to 
Standard Normal Distribution for Angles of Incidence at 
25.8 Jet Diameters Downwind from the Nozzle 
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Comparison of Horizontal Distribution of Temperature to 
Standard Normal Distribution for Angles of Incidence at 
51.7 Jet Diameters from the Nozzle 
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Comparison of Horizontal Distribution of Temperature to 
Standard Normal Distribution for Angles of Incidence at 
86.2 Jet Diameters Downwind from the Nozzle 
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Comparison of Horizontal Distribution of Tempera 
Standard Normal Distribution for Angles of Inci 
120.6 Jet Diameters Downwind from the 
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20 30 Downwind Distance (JD) 70 80 

Figure 33 

Spatial Representation of Plume Width as a Function of 
Downwind and Horizontal Distances from the eid of the 
Nozzle for Different Angles of Incidence to the Wind 
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Figure 34 

Distance Downwind and Normal to Wind from the Nozzle 
Where the Plume is Uniformly Distributed 
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Horizontal Plume Width vs. Downwind Distance 



Wind on (about 40 ft/sec) 



xj 

d 

4-4 

^ 3 
0 rO 




O 

eg 



o 

o 



✓— \ 
d 

o 

c 

d 

4-4 

CO 



Q 

xJ 

d 

’rl 

£ 



O 

Q 



vO 

CO 

d 

d 

&0 

•H 

Pli 



O 



o 

CN 



d 

a 

d 

d 

4J 

CO 

• 1-4 

Q 



X) 

C 



s 



co 

> 



x: 

4J 

X) 

•H 



d 

a 



V4 

d 

> 



62 



700 

550 



CO CO 
CD Q) 

Q> <D 
U V-i 
DO DO 
0) 0) 
O Q 



II II 



CD 0) 
M H 
3 d 

4J 4J 

cd cd 

o) a) 
Pu Ok 

e e 

a) <d 
H H 

» I. 



I i 



— P 

co 





(*UT) eOUE^STQ JH0T^J9A 



o 

CM 



O 

o 



o 

G\ 



O 

00 



o 

r^> 



O ^ 
VO Q 



O 

LO 



o 




o 

CO 



o 

CM 



* 



63 



Plume Rise vs. Downwind Distance 
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Jet Penetration vs. Downwind and Horizontal Distance from 
the Nozzle for Different Angles of Incidence to the Wind 




Jet Axis Penetration Distance (JD) 
Figure 39 



Jet Penetration Distance along the Jet 
Axis vs. Angle of Incidence 
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Figure 40 

1-Sigma Width of Plume At Gaussian Position 
vs. Angle of Incidence to the Wind 
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